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We report measurements of branching fractions and charge asymmetries of exclusive decays of 
neutral and charged B mesons into two-body final states containing a charmonium state and a light 
strange meson. The charmonium mesons considered are J/tp , tp{2S) and Xci, and the light meson is 
either K or K* . We use a sample of about 124 million BB pairs collected with the BABAR detector 
at the PEP-II storage ring at the Stanford Linear Accelerator Center. 

PACS numbers: 13.25.Hw, 12.15.Hh, 11.30.Er 



Nonleptonic decays of B mesons provide tests of both 
strong- and weak-interaction dynamics. Decays B — > 
{cc)K^*^ are particularly illuminating as they involve 
three kinds of mesons: one with a heavy quark and a 
hght quark, one with two heavy quarks, and one with two 
light quarks. Phenomenological models (see Refs.[2-12] 
of Ref. 0) give estimates for branching fractions and for 
ratios of the decays to K and K*. Some branching frac- 
tions and ratios have been reported but others have not, 
so more stringent tests of the models are possible. The 
Standard Model predicts small differences between the 
branching fractions for positive and negative B mesons, 
i.e. small direct CP violation Large charge asymme- 
tries would be evidence for new physics. Limits on direct 
CP violation would constrain extensions of the Standard 
Model. Very few measurements have been reported in 
B {cc)K^*^ modes 0. The decay processes studied in 
this paper are listed in Tabled 

TABLE L Decay processes studied in this analysis. The (cc) 
resonance is either a J/tp , ip{2S), or xa- For all processes, 
when relevant, we use the secondary decay modes J/tp ,tp{2S) 
£+£-, Xci J/ip 7, and K°s -> tt+tt", tv° ^ 77. 



Decay Channel 


Secondary K* decay mode 


B" (cc) K*" 
B+ {cc) K'+ 
B° (cc) K'^ 
B+ (cc) K+ 


A-u IT-, Ki 
K'+ ^ K+ n°, Kl n+ 



The data sample used in this analysis contains 124 mil- 
lion BB events collected with the BABAR detector at the 
PEP-II asymmetric e+e^ storage ring. This represents 
a total integrated luminosity of 112.4 fb^^ taken on the 
T(4S) resonance. The BABAR detector is described in de- 
tail elsewhere Surrounding the interaction point, a 
five-layer double-sided silicon vertex tracker (SVT) pro- 
vides precise reconstruction of track angles and B de- 
cay vertices. A 40-layer drift chamber (DCH) provides 
measurements of the transverse momenta of charged par- 
ticles. An internally reflecting ring-imaging Cherenkov 
detector (DIRC) is used for particle identification. A 
CsI(Tl) crystal electromagnetic calorimeter (EMC) is 
used to detect photons and electrons. The calorimeter is 
surrounded by a 1.5-T magnetic field. The flux return is 
instrumented with resistive plate chambers (RPC) used 
for muon and neutral-hadron identification. 



Multihadron events are selected by demanding a mini- 
mum of three reconstructed charged tracks in the polar- 
angle range 0.41 < < 2.54 rad, where is defined 
in the laboratory frame. Charged tracks must be re- 
constructed in the DCH and are required to originate 
within 1.5 cm of the beam in the plane transverse to it 
and within 10 cm of the beamspot along the beam di- 
rection. Events are required to have a primary vertex 
within 0.5 cm of the average position of the interaction 
point in the plane transverse to the beamline, and within 
6 cm longitudinally. Charged tracks are required to in- 
clude at least 12 DCH hits and to have a transverse mo- 
mentum pt >100 MeV/c. Photons are reconstructed 
from EMC clusters. The lateral energy profile (LAT) 
is used to discriminate electromagnetic from hadronic 
clusters. Photons are required to have a minimum en- 
ergy of 30 MeV, to satisfy LAT < 0.8, and to be in the 
fiducial volume 0.41 < 9 < 2.41 rad. Electron candi- 
dates are selected using information from the EMC (LAT 
and Zernike moment A42 0), the ratio of the energy 
measured in the EMC to the momentum measured by 
the tracking system, the energy loss in the drift cham- 
ber, and the Cherenkov angle measured in the DIRC. 
Electrons are also required to be in the fiducial volume 
0.41 < 9 < 2.41 rad. Muon candidates are selected using 
information from the EMC (energy deposition consistent 
with a minimum ionizing particle) and the distribution of 
hits in the RPC. Muons are required to be in the fiducial 
volume 0.3 < 9 < 2.7 rad. We select charged kaon and 
pion candidates using information from the energy loss in 
the SVT and DCH, and the Cherenkov angle measured 
in the DIRC. Kaon candidates are required to be in the 
fiducial volume 0.45 < 9 < 2.45 rad. 

The selection has been optimized by maximizing the 
ratio S/\/S + B, where S and B are the number of 
expected signal and background events obtained from 
Monte Carlo simulation. The J/tp candidates are re- 
quired to have an invariant mass 2.95 < Mg+g- < 3.14 
GeV/c2 or 3.06 < M^+^- < 3.14 GeV/c^ for J/i/j 
e+e^ or J/tp M^M^ decays respectively. The i'{2S) 
candidates are required to have an invariant mass 3.44 < 
Afg+g_ < 3.74 GeV/c^ or 3.64 < M^+^- < 3.74 GeV/c^ 
for ip{2S) — > e+e^ or ip{2S) A^^M^ decays respec- 
tively. For J/^p,ip{2S) e+e^ decays, electron candi- 
dates are combined with photon candidates in order to 
recover some of the energy lost through brcmsstrahlung. 
In the Xci reconstruction, J/ip candidates are selected as 
described above. The associated 7 has to satisfy LAT < 
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state. To determine the extent of peaking background, 
the niES distribution for simulated BB events is fitted 
within the AE signal region by an ARGUS function and 
a Gaussian function. The peaking background is taken 
as the integral of the Gaussian portion. 
The branching fractions are obtained as: 



BF : 



Ns 



N 



BB 



X e X / 



(1) 



0.8, A42 <0.15 and has to have an energy greater than 
0.15 GeV. The Xci candidates are required to satisfy 
0.35 < M(>+e-^ - M(>+f>- < 0.45 GeV/c^, where £ repre- 
sents an electron or a muon. The tt*^ —> 77 candidates 
are required to satisfy 0.113 < M^^ < 0.153 GeV/c^. 
Both photons have to satisfy LAT < 0.8. The energy of 
the soft photon has to be greater than 0.050 GeV and 
the energy of the hard photon has to be greater than 
0.150 GeV. The K'^ 7r+ vr^ candidates are required 
to satisfy 0.489 < M^+^- < 0.507 GeV/c^. In addition, 

the K° flight distance defined as the distance between where N^j^ is the number of BB events, e is the selection 
the reconstructed B and K'^ vertices must exceed 1 mm, 
and the angle between the momentum and its flight 
direction in the plane transverse to the beam axis must 
be less than 0.2 rad. The K*° and K*+ candidates are 
required to satisfy 0.796 < Mk^ < 0.996 GeV/c^ and 
0.792 < MKn < 0.992 GeV/c^, respectively. In addition, 
for the sake of suppressing background from events with 
soft pions, for channels having a 7r° in the final state, the 
cosine of the angle between the K momentum and the B 
momentum in the K* rest frame has to be less than 0.8. 

The B candidates are reconstructed by combining 
charmonium and kaon candidates and are characterized 
by two kinematic variables: the difference between the 
reconstructed energy of the B candidate and the beam 
energy in the center-of-mass frame AE = E^ — El^^^, 
and the beam energy-substituted mass rriESi defined as 
mES = V ^beam ^ P*B ' where the * refers to quantities 



in the center-of-mass and ps is the B momentum. For a 
correctly reconstructed B meson, AE is expected to peak 
at zero and the energy-substituted mass mES at the B 
meson mass, 5.279 GeV/c^. Only one reconstructed B 
meson is allowed per event. For events that have multi- 
ple candidates, the candidate having the smallest \AE\ 
is chosen. The analysis is performed in a region of the 
Toes vs AE plane defined by 5.2 < toes < 5.3 GeV/c^ 
and -0.12 < AE < 0.12 GeV. A AE channel-dependent 
signal region is subsequently defined. The toes distribu- 
tions within the AE signal region for candidate events 
are shown on Figure ^ 

There are two components to the residual background 
in the AE signal region: the combinatorial background 
and a peaking component (component of the background 
that peaks at the same values of AE and toes as the sig- 
nal). The number of signal events Ns is determined from 
the number of candidate events, iVcandi after subtracting 
the peaking background. For this purpose, the toes dis- 
tribution within the AE signal region is fitted to the sum 
of an ARGUS function |^ , which models the combinato- 
rial background, and a Gaussian function. The value of 
A^cand is givcn by the integral of the Gaussian component. 
There are two contributions to the peaking background. 
The first is the cross-feed component that is due to B 
(cc) K* events from the three channels other than 
the one under consideration. The second contribution is 
from other B decays with a J/tp or ■0(25') hi the final 



BB 

efficiency and / is the total secondary branching fraction. 
For channels with a K* in the final state, the cross-feed 
contribution depends on the branching fractions that are 
being measured. It was estimated by an iterative proce- 
dure and found to be small. 

The systematic errors arise from the uncertainty on 
the number of BB events (1.1%), the secondary branch- 
ing fractions (taken from Ref. 8]), the estimate of 
the selection efficiency, and the knowledge of the back- 
ground. For the tracking efficiency, an error of 1.3% per 
track has been used. For the particle identification ef- 
ficiency, the systematic error varies between 0.2% and 
3.7%. The uncertainty on the detection and energy mea- 
surement of photons is 2.5%, common to all channels, 
plus a small channel-dependent correction. The uncer- 
tainty on the 7r° reconstruction is 5.0% for all chan- 
nels, plus a channel-dependent correction. The over- 
all selection efficiency depends on the angular distribu- 
tion. The efficiency can be written as e = a + Aob, 
where a and b are obtained from the K* helicity an- 
gle by a = 3/4 /(I - cos'^ 0K')e{0K') sin{0K')deK'' and 
6 3/4 /(3cos2 Ok' - l)e{eK')sm{eK')deK', and Aq is 
the fraction of the longitudinal K* polarization The 
values of a and b are obtained from simulation. A sys- 
tematic error varying from 3.4% to 8.6% has been derived 
for the branching fractions. In the default fit, the shape 
parameter of the ARGUS function is not constrained. 
To determine a systematic error due to the combinato- 
rial background, a second fit with the shape parameter of 
the ARGUS function fixed to the value obtained from fit- 
ting the data in the AE sideband region was performed. 
The systematic uncertainty on the combinatorial back- 
ground has been taken as 50% of the difference between 
the number of events obtained from the default fit and 
from the second fit. For the cross-feed component to the 
peaking background, the uncertainty of the correspond- 
ing branching fractions, taken from Ref. 0, has been 
assigned as the systematic error. For the contribution 
coming from other B decays with a J/ip or a ip{2S) in 
the final state, a 50% error has been assigned, accounting 
for the poor knowledge of the branching fractions of the 
contributing decay modes. Overall, the dominant con- 
tribution to the systematic error is from the secondary- 
decay branching fractions in the case of i'{2S) K'^*^ and 
Xci K, the efficiency determination in the case of J/ip 
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FIG. 1: mES distributions and fits within the AE signal 
region. From left to right, the columns show the distributions 
for the J/i> , ip{2S) and Xci channels. From top to bottom, 
the rows show the distributions for the B — > (cc) if*" {K^ 
7r°), B (cc) K*° {K+ TT"), B ^ (cc) K*+ {Kl 7r+), B ^ 
(cc) K*+ {K+ 7r°), B -> (cc) Tfg and B (cc) A'+ decay 
modes. The dashed and dotted lines show the combinatorial 
and peaking backgrounds, respectively. 



TABLE II: Measured branching fractions for exclusive decays 
of B mesons to charmonium and kaon final states. The first 
error is statistical and the second systematic. 



Channel 



Branching fraction ( x 10 



B" 




13.09±0.26±0.77 




J/ip K 


i4.54±U.47±U.97 


B+ 


-> J/iP K+ 


10.61±0.15±0.48 


B° 


J/iP K° 


8.69±0.22±0.30 


B" 


iP{2S) K"" 


6.49±0.59±0.97 


B+ 


tP{2S) K*+ 


5.92±0.85±0.89 


B+ 


iP{2S) K+ 


6.17±0.32±0.44 


B° 


Tp{2S) K° 


6.46±0.65±0.51 


B° 


- Xci K'' 


3.27±0.42±0.64 


B+ 


- Xci K*+ 


2.94±0.95±0.98 


B+ 


^ Xci K+ 


5.79±0.26±0.65 


B" 


Xci K 


4.53±0.41±0.51 



can compute R^/'^. Using the ratio of the charged to 
neutral B meson lifetimes t^+Itbo = 1.086 ± 0.017 
we obtain: 



T{r{AS) B+B-) 

r(r(45') ^ SOB") 



1.06 ±0.02 ±0.03, (2) 



where the first error is statistical and the second sys- 
tematic. The ratios of the branching fractions for B — > 
{cc)K* and B {cc)K for the three (cc) states are pre- 
sented in Table UTTI For each of the charmonium states, 
the average of the charged and neutral measurements is 
also shown. 



TABLE III: Results for ratios of the branching fractions for 
charmonium and a K* versus charmonium and a K. The first 
error is statistical and the second systematic. 



Ratio 




Result 


B{B" ^ J/ipK'"')/B{B" ^ J/^K") 
B{B+ .J/iPK''+)/B{B+ J/ipK+) 
B{B -> J/iPK*)/B{B J/iPK) 


1.51±0.05±0.08 
1.37±0.05±0.08 
1.44±0.04±0.06 


B{B" -> iP(2S)K*")/B{B" ^ '>pi2S)K") 
B{B+ i>{2S)K'-+)/BiB+ ^{2S)K+) 
B{B ^ tP{2S)K*)/B(B ip{2S)K) 


1.00±0.14±0.09 
0.96±0.15±0.09 
0.98±0.10±0.07 


B{B" 
B{B+ 
B{B - 


^ XciK'^)/B{B" -> xciif") 
^ XciK*+)/B{B+ XciK+) 
-> XciK')/B{B ^ xcii^) 


0.72±0.11±0.12 
0.51±0.17±0.16 
0.65±0.09±0.10 



K'^*\ and the background subtraction in the case of Xci 
K*. 

The branching fraction measurements are summarized 
in Table CH Wc have taken the ratio T{T{AS) 

B+B-)/T{r{4:S) B^B°), to be unity for the deter- 
mination of the branching fractions. Assuming isospin 
invariance in the B — > charmonium K (K*) decays, we 



Finally, we have measured the charge asymmetries 

^ B{B+ ^ (cc)j^+(*)) - B{B- ^ {cc)K~(*^) 

~ B{B+ ^ (cc)i^+(*)) + B{B- ^ (cc)if -(*)) ' ^ ' 

using efficiencies determined separately for the two 
charges. The results are presented in Table Hvl No sta- 
tistically significant asymmetry is observed. 
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TABLE IV: Results for charge asymmetries. The first error 
is statistical and the second systematic. 



Final state 



Asymmetry 



J/V- K*+ 
iP{2S) K+ 
ip[2S) K*+ 
Xci K+ 



■0.030±0.014±0.010 
0.048±0.029±0.016 
0.052±0.059±0.020 
■0.077±0.207±0.051 
0.003±0.076±0.017 
■0.47I±0.378±0.268 
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In summary, branching fraction measurements of ex- 
clusive B decays to charmonium (J/'0, ^(25*) and Xci) 
and K or K* have been presented. Our results for J/tp 
and 'ip{2S) are in good agreement with previous measure- 
ments !^ and exhibit comparable or superior precision. 
Our Xci results have much better precision than earlier 
measurements. The B'^ —> Xci K*'^ mode was previ- 
ously unmeasured. Assuming isospin invariance, we find 
the ratio of charged- to neutral- i? meson production on 
the T(4S) resonance to be compatible with unity within 
1.7 standard deviations. No direct CP violation has been 
observed in the charge asymmetries. 
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